Controllability of the bilinear Schrodinger equation with several 
controls and application to a 3D molecule* 
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Abstract — We show the approximate rotational controllabil- 
ity of a polar linear molecule by means of three nonresonant 
linear polarized laser fields. The result is based on a general 
approximate controllability result for the bilinear Schrodinger 
equation, with wavefunction varying in the unit sphere of 
an infinite-dimensional Hilbert space and with several control 
potentials, under the assumption that the internal Hamiltonian 
has discrete spectrum. 



I. Introduction 

Rotational molecular dynamics is one of the most im- 
portant examples of quantum systems with an infinite- 
dimensional Hilbert space and a discrete spectrum. Molec- 
ular orientation and alignment are well-established topics 
in the quantum control of molecular dynamics both from 
the experimental and theoretical points of view (see [19], 
[20] and references therein). For linear molecules driven by 
linearly polarized laser fields in gas phase, alignment means 
an increased probability direction along the polarization 
axis whereas orientation requires in addition the same (or 
opposite) direction as the polarization vector. Such controls 
have a variety of applications extending from chemical reac- 
tion dynamics to surface processing, catalysis and nanoscale 
design. A large amount of numerical simulations have been 
done in this domain but the mathematical part is not yet 
fully understood. From this perspective, the controllability 
problem is a necessary step towards comprehension. 

We focus in this paper on the control by laser fields of 
the rotation of a rigid linear molecule in K 3 . This con- 
trol problem corresponds to the control of the Schrodinger 
equation on the unit sphere S 2 . We show that the system 
driven by three fields along the three axes is approximately 
controllable for arbitrarily small controls. This means, in 
particular, that there exist control strategies which bring 
the initial state arbitrarily close to states maximizing the 
molecular orientation [21]. 
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A. The model 

We consider a polar linear molecule in its ground vi- 
bronic state subject to three nonresonant (with respect to 
the vibronic frequencies) linearly polarized laser fields. The 
control is given by the electric fields E = {u\, 1*2,^3) 
depending on time and constant in space. We neglect in 
this model the polarizability tensor term which corresponds 
to the field-induced dipole moment. This approximation is 
correct if the intensity of the laser field is sufficiently weak. 
Despite its simplicity, this equation reproduces very well 
the experimental data on the rotational dynamics of rigid 
molecules (see [20]). 

Up to normalization of physical constants (in particular, in 
units such that H = 1), the dynamics is ruled by the equation 



at 



= — Atp(9, ip, t) + (ui (i) sin cos ip 

U2(t) sin 8 sin <p + Us(t) cos9)ip(6, <p, t) 

(1) 



where 6, <p are the spherical coordinates, which are related 
to the Euclidean coordinates by the identities 



sinf^cos^, y = sinf? sin 95, 



cos 9, 



while A is the Laplace-Beltrami operator on the sphere 
(called in this context the angular momentum operator), i.e., 



A 



1 d 



sin t 



sin 9 



d 



+ 



1 



d 2 



The wavefunction 

H = L 2 (S 2 ,C). 



09 J ' sin 2 9 dip 2 ' 
t) evolves in the unit sphere S of 



B. The main results 

In the following we denote by ip(T; ipo, u) the solution at 
time T of equation (HJ, corresponding to control u and with 
initial condition ip(0; ipo, u) — tpo, belonging to S. 

Our main result says that (Q3 is approximately controllable 
with arbitrarily small controls. 

Theorem 1.1: For every ip°, ip 1 belonging to >S and 
every e, 61,62,63 > 0, there exist T > and u 6 
^([OjTjjrVi] x [0,5 2 ] x [0,63]) such that H^ 1 - 

V>Cr ; VV)|| < £ . 

The proof of the result is based on arguments inspired by 
those developed in [11], [7]. There are two main difficulties 
preventing us to apply those results to the case under consid- 
eration: firstly, we deal here with several control parameters, 
while those general results were specifically conceived for 
the single-input case. Notice that, because of symmetry 
obstructions, equation (fl]i is not controllable with only one 



of the three controls u%, 1*2, 113. Secondly, the general theory 
developed in [11], [7] is based on nonresonance conditions 
on the spectrum of the drift Schrodinger operator (the inter- 
nal Hamiltonian). The Laplace-Belatrami operator on S 2 , 
however, has a severely degenerate spectrum. It is known, 
indeed, that the l-\h eigenvalue — i£(£ + 1) has multiplicity 
2£ + 1. In [11] we proposed a perturbation technique in 
order to overcome resonance relations in the spectrum of 
the drift. This technique was applied in [8] to the case of the 
orientation of a molecule confined in a plane driven by one 
control. The planar case is already technically challenging 
and a generalization to the case of three controls in the space 
will hardly provide an apophantic proof of the approximate 
controllability result. We therefore provide a general multi- 
input result which can be applied to the control problem 
defined in (|T), up to the computation of certain Lie algebras 
associated with its Galerkin approximations. 

The structure of the paper is the following: in the next sec- 
tion we present the general multi-input abstract framework 
and we recall some previously known controllability and 
non-controllability results. In Section [III] we prove our main 
sufficient condition for approximate controllability. Finally, 
in Section [IV] we prove that the abstract result applies to 
system ((TJ. 

II. Abstract framework 

Definition 2.1: Let T-L be an infinite-dimensional Hilbert 
space with scalar product (•, •) and A, B%, . . . ,B p be (pos- 
sibly unbounded) linear operators on "H, with domains 
D(A),D(B\), D(B p ). Let U be a subset of W. Let us 
introduce the controlled equation 



~dt 



(t) = (A+u 1 {t)B 1 + - ■ ■+u p (t)B p )ijj(t),u(t) G U C 



(2) 

We say that (A, B\, . . . , B p , U, <£) satisfies (21) if the fol- 
lowing assumptions are verified: 

(211) $ = (<frk)keN is an Hilbert basis of H made of 
eigenvectors of A associated with the family of 
eigenvalues (iAfc) feeN ; 

(212) cf) k G D(Bj) for every k G N, j = 1, . . . ,p; 

(213) A + u 1 B 1 + h u p B p : span-j> & | k G N} -> U 

is essentially skew-adjoint for every u G U; 

(214) if j ^ k and Aj = A fc then {fa , B^) = for 
every I — 1 , . . . , p. 

If (A, Bi,...,B p , U, 3?) satisfies (21) then, for every 
(til, . . . , Up) G M p , A + u\B\ + ■ ■ ■ + UpBp generates a 



unitary group e 



t(A+u 1 B 1 + ---+u p B v ) 



It is therefore possible 



to define the propagator at time T of system (Q~|) 
associated with a p-uple of piecewise constant controls 
u(t) = (ui(t), . . . ,u p (t)) by concatenation. If, moreover, 
the potentials B\ , . . . , B p are bounded operators then the 
definition can be extended by continuity to every L°° control 
law. 

Definition 2.2: Let (A, B, U, $) satisfy (21). We say that 
(O is approximately controllable if for every i/jq, in the 
unit sphere of H and every e > there exist a piecewise 



constant control function u : [O.T] — > U such that \\1p1 — 
T«(Vo)||<£. 

Definition 2.3: Let (A, B, U, $) satisfy (21) . We say that 
(ffj is approximately simultaneously controllable if for every 
r in N, ipi, . . . , t/v in %, T in U(H), and e > there exists 
a piecewise constant control u : [0,T] — > U such that 



< 



1, 



A. Short review of controllability results 

The controllability of system (J2j is a well-established topic 
when the state space H is finite-dimensional (see for instance 
[12] and reference therein), thanks to general controllability 
methods for left-invariant control systems on compact Lie 
groups ([10], [14]). 

When H is infinite-dimensional, it is known that the bilin- 
ear Schrodinger equation is not controllable (see [2], [22]). 
Hence, one has to look for weaker controllability properties 
as, for instance, approximate controllability or controllability 
between eigenstates of the Schodinger operator (which are 
the most relevant physical states). In certain cases where 
the dimension of the domain where the controlled PDE is 
defined is equal to one a description of the reachable set 
has been provided [3], [4], [5]. For dimension larger than 
one or for more general situations, the exact description 
of the reachable set appears to be more difficult and at 
the moment only approximate controllability results are 
available. Most of them are for the single-input case (see, 
in particular, [6], [7], [11], [15], [16], [18], [17]), except for 
some approximate controllability result for specific systems 
([13]) and some general approximate controllability result 
between eigenfunctions based on adiabatic methods [9]. 

B. Notation 



Set b 



(0 
> 



4>j,Bi<pk), I = 1, • ■ • ,p- For every n in N, 
define the orthogonal projection 

j<n 

Given a linear operator Q on TL we identify the linear 
operator n n Q^n preserving span{0i, . . . , </>„} with its n x 
n complex matrix representation with respect to the basis 

((f>l,...,<f> n ). 

III. Main abstract controllability result in the 

MULTI-INPUT CASE 

Let us introduce the set of spectral gaps associated 
with the TV-dimensional Galerkin approximation as 

Ejv = - Afc I I j,k = 1,. . . , N,\j ^ A fc }. 

For every a G Sjv. let 

bW( Vi , . . .,v p ) jik = Mf»+. . .^U,iva»|. 

The N x N matrix B„'(vi,...,v p ) corresponds to the 
choice of the controls v\ , . . . , v p and to the "activation" of 
the spectral gap a. Define 

M N = {Bi N \ Vl ,...,v p ) I <7G E N , Vl ,...,v p G [0,1]} 



and 
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M e su(n) ViV >n3Q e M N s.t. Q 

The set A^q represents "compatible dynamics" for the n- 
dimensional Galerkin approximation (compatible, that is, 
with higher dimensional Galerkin approximations). 

Theorem 3.1 (Abstract multi-input controllability result): 
Let U = [0,S] P for some S > 0. If for every n G N there 
exist n > rig such that 

LieMft =Bu(ra), (3) 

then the system 

x = (A + Ui-Bi + • • • + u p Bp)x, u G [/, 
is approximately simultaneously controllable. 

A. Preliminaries 

The following technical result, which we shall use in the 
proof of Theorem 13.11 has been proved in [7]. 

Lemma 3.2: Let k be a positive integer and 71, . . . , 7* G 
R \ {0} be such that -yi | ^ |^ 1 for j = 2, . . . , k. Let 

Vj(t) = (e«^,... J e < *T-)- 
Then, for every r G R, we have 



convv?([T ,oo)) D i/S 1 x {(0,...,0)}, 



where = 11^2 cos (sfc) > ^- Moreover, for every R > 
and £ G S 1 there exists a sequence (tk)k£fi such that 
tk+i -tk> R and 



1 ' l 

, lim T/J = ( v &°> 



,0). 



fc=i 



B. 7Ime reparametrization 

For every piecewise constant function z(t) 



Y,k=i z kX[s k _ 1 , Sk ){t) such that z k > 0, for every 

k = 1,..., if, and ^(t) = Efc s i"fc' ) X[a»_ 1 ,. fc )(*) 
with j = 1, . . . ,p, we consider the system 



-^0) = (2(*)A + Vl (i)Bi + 



■v p (t)B p )^(t). (4) 



System (0]) can be seen as a time-reparametrisation of 
system (f2]). Let ^(t) be the solution of (f2]i with initial condi- 
tion tpQ G H associated with the piecewise constant control 



u(-) with components m^O = EfcLi u fc J X[t k - lt t„)(-)> 3 = 
1, . . . ,p. If s fc = + s k -u s = 0, vjp = u ( k j) z k for 

every k = 1, . . . ,K, j = 1, . . . ,p, then the solution ip(t) 
of © with the initial condition i/jq G % associated with the 
controls • • • , satisfies 



Jo^ X[tk - Utk)(S)dS ) 



Controllability issues for system <J2J and (@J are equivalent. 
Indeed, consider piecewise constant controls z : [0,^] — >• 

[1/(5, OO), z(i) = EfcLl^fcX[a fc _ 1)afc )(*) and V 3 : [°' T "] ~* 

[0,1], Vj (t) = E£=i«£°X[.*-i..*)(*) with j =l,...,p, 
achieving controllability (steering system (|4]i from ipj to 
Ttpj, j = 1, . . . , r in a time T„) . Then the controls = 



,0'), 



,0') 



EfcLi "* ^[t*-!,**)' i = • • - >P defined by u { k o) = v^/z k 
and to = Q,t k = (fife — Sk-i)Zk+tk-i, steer system (O from 
ipj to Ti/>j, i = 1, • • • j T in a time T u . 

C. Interaction framework 

Let u(t) — f Q z(s)ds, and Wj(t) = J Q t Vj(s)ds for j = 
I,..., p. Let ij)(f) be the solution of (HJi with initial condition 
ipo G H associated with the controls z(t), i>i(t), . . . , v p (t) 
and set 

y(t) = e-^ A m- 



For ui, vt, . . . , v p G R set 0(w, vi, . . . , v p ) 
■ ■ ■ + Vp(t)B p )e u ' A , then y(t) satisfies 



y(t) = e(cj(t),v 1 (t),...,v p (t))y(t). 



(5) 



Note that 

0(w, 



Vpfjk 



k, Q(u>,Vi, . ■ -,V p )cpj) 



+ w D 6 



(p) 



Notice that |y(i)| = \ip(t)\, for every t G [0,T V ] and 
for every (p + l)-uple of piecewise constant controls z : 

[0,T V ] -> [1/(5, +00), ui,..., up : [0,T„] ^ [0,1]. 

Z3. Galerkin approximation 

Definition 3.3: Let TV G N. The Galerkin approximation 
of (0 of order TV is the system in H 



(6) 



where 6 (A,) (a;, Ui, . . . , v p ) = n N Q(u,Vi, . . . , v p )tt n 



(S(uj,v 1 ,...,v p ) jk )^ k=1 . 



E. First step: choice of the order of the Galerkin approxi- 
mation 

In order to prove approximate simultaneous controllability, 
we should take r in N, ipi, . . . , ip r in %, T in XJ(H), and 
e > and prove the existence of a piecewise constant control 
u : [0, T] -)• £/ such that 



< e, 



L 



Notice that for no large enough there exists U G SU(no) 
such that 

\{(j)j,T^ k ) - (7T no (^-,Z77r no V>fc)| < e 

for every 1 < A; < r and j G N. This simple fact suggest to 
prove approximate simultaneous controllability by studying 
the controllability of © in the Lie group SU{no). 



F. Second step: control in SU(n) 

Let n > no satisfy hypothesis (01. It follows from standard 
controllability results on compact Lie groups (see [14]) that 
for every U G SU(n) there exists a path M : [0, T v ] -> Mq 
such that 



exp / M(s) ds = U, 
Jo 

where the chronological notation exp J Q V s ds is used for the 
flow from time to t of the time-varying equation q = V s (q) 
(see [1]). More precisely, there exists a finite partition in 
intervals of [0,T„] such that for every t G Ik either 

there exist v\, . . . , v p G [0, 1] and a 6 such that 

M{t) =7r„flW » p )7r n , 

or 

M(t)=A(")-^/„. 

n 

In particular, 

M(% fc = 0, for every i G [0, T„], j <n,k> n. (7) 

G. Third step: control of A4 jv 

Lemma 3.4: For every iV G N, £ > 0, and for ev- 
ery piecewise constant ^i,...,^ : [0,T„] — > [0,1] and 
a : [0,T„] — > there exists a sequence (^/i(-))h,eN of 
piecewise constant functions from [0,T V ] to [1/5, oo), such 
that 



<d (N \z h (s),vi, . . . ,v p )ds 



-> 



uniformly with respect to t G [0, T v ] as h tends to infinity. 

In other words, every piecewise constant path in M. n can 
be approximately tracked by system ©. 
Proof. Fix N G N. We are going to construct the control 
Zh by applying recursively Lemma 13.21 Consider an interval 
[tk,tk+i) in which Vj(t), j — l,...,p, and a{t) are 
constantly equal to Vj G [0, 1], j = 1, . . . ,p, and a G Sat 
respectively. Apply Lemma l3~2l with 71 = a, {72, . . . , 7 K } = 
Sat\{it}, R~T and To = To(fe) to be fixed later depending 
on k. Then, for every 77 > 0, there exist h = h(k) > 1/rj and 
a sequence (iu„)„ =1 such that w\ > to, — w 1 ^ 
and such that 



> R, 



h 

IE' 

Q = 1 



,'i(A z -A m )i 



-(TV) 



(iV) 



J <T,\Xl-\r, 



< 7 h 



Set t£ = t fc + (t fc+1 - t k )a/h, a 
the piecewise constant function 

h{k) 



0, . . . , h, and define 



(8) 



Note that by choosing ro(fc) = + R for fc > 1 and 

T o(0) = R we have that Wfc(i) is non-decreasing. 

Following the smoothing procedure of [7, Proposition 5.5] 
one can construct the desired sequence of control z^. The 
idea is to approximate U!h(t) by suitable piecewise linear 
functions with slope greater than 1/5. Then Zh can be 
constructed from the derivatives of these functions. □ 

As a consequence of last proposition by [1, Lemma 8.2] 
we have that 



exp 



e^(z h (s),v 1 (s),...,v p (s)ds 



exp / S^ ) ) (wi(s),...,u p (s)) 
Jo 



ds 



uniformly with respect to t G [0, T v ] as h tends to infinity. 

H. Fourth step: control of the infinite-dimensional system 

Next proposition states that, roughly speaking, we can 
pass to the limit as N tends to infinity without losing 
the controllability property proved for the finite-dimensional 
case. Its proof can be found in [7, Proposition 5.6]. It is 
based on the particular form (fT) of the operators involved, 
since the fact that the operator has several zero elements 
guarantees that the difference between the dynamics of the 
infinite-dimensional system and the dynamics of the Galerkin 
approximations is small. 

Proposition 3.5: For every e > 0, for every 5 > 0, and for 
every trajectory U G SU(n) there exist piecewise constant 
controls Uj : [0,T U ] — > [0,5},j — l,...,p such that the 
associated propagator T u of © satisfies 



(TTn^i ) UTr n (f)) I — I ((pj , Ty u cj>) I 



< £ 



J with 



= 1 and for every 



for every (j) G span{0i, . . 
j in N. 

We recall now a controllability result for the phases 
(see [7, Proposition 6.1 and Remark 6.3]). This property, 
stated in the proposition below, together with the controlla- 
bility up to phases proved in the previous section, is sufficient 
to conclude the proof of Theorem 13.11 

Proposition 3.6: Assume that, for every T G XJ(H), m 
in N, 5 > 0, and e > 0, there exist T u > and piecewise 
constant controls Uj : [0, T u ] — > [0, 5], j = 1, . . . ,p such that 
the associated propagator T u of equation (f2]i satisfies 



m,t<i>)\-m,nj)\ 



< e, 



k>0 a = l 



for every j G N and tf> G span{0i, . . . , </> m } with \\<f>\\ = 1. 
Then (fJJ is simultaneously approximately controllable. 

IV. 3D MOLECULE 

Let us go back to the system presented in the introduction 
for the orientation of a linear molecule, 

ihip = —Aip + (u\ cos 9 + U2 cos^sin^ + U3 sin ip sin 6)ip, 

(9) 

where ^)e« = L 2 (S 2 ,C). 



A basis of eigenvectors of the Laplace-Beltrami operator 
A is given by the spherical harmonics YJ n (6,(p), which 
sastisfy 

AY e m (6,<p) = -t(£ + l)Y e m (e,<p). 

We are first going to prove that for every £ S N the system 
projected on the (A£ + 4)-dimensional linear space 



C := spanjY^, 



,Yp , Y» 



Y e+u 
■ ) I t+i i 



UeM 



41+4 




A. Matrix representations 

Denote by Jn the set of integer pairs {(j, k) j — £, £ + 
1, k — —j, . . . , j}. Consider an ordering u) : {1, . . . , A£ + 
4} — > Jg. Let ej^k be the (4^+4)-square matrix whose entries 
are all zero, but the one at line j and column k which is equal 
to 1. Define 



By a slight abuse of language, 
e^fc. The analogous identification can be used to define 



also set e u(i ) iU)(fe ) 



E, 



Thanks to this notation we can 
the matrices corresponding to the 
(projected on C). A computation 
potential in the z direction, —i cos 
a matrix representation with respect 



conveniently represent 
controlled vector field 
shows that the control 
9, projected on £, has 
to the chosen basis 



B3 — ^ P(.,m,E(t,m),{i+l,m) 
m=-l 



with 



P(.: 



! (l + l) 2 -m 2 

(2e + i)(2e + 3)' 



Similarly, we associate with the control potentials in the 
x and y directions, — icosipsmO and — ismipsmO respec- 
tively, the matrix representations 

I 



B2 — ^ (<tt,m^(*,m),(&fl,m-l 



) + Qe,—mE(t,m),{e+l,m+l)) 



where 



(£-m + 2)(£-m + l) 
4(2£ + 1)(2£ + 3) ' 



The matrix representation of the Schrodinger operator iA 
is the diagonal matrix 



where 



' 1 e+i 

is controllable. More precisely, chosen a reordering ((j)k)keN 
of the spherical harmonics in such a way that 

{<f>k\k=l,..., 4^+4} = {Yf l , Yf, Y e -\-\ . . . , 

we are going to prove that 



*(j,k) = + 1), for (j, k) e Je. 



Now consider A = A — I, in such a way that 



4(£+l) J 

tx{A) = 0. Hence, A = Y,(j, k )eJ e a m e m,{o,k) 
2£ + 3 



where 



-, for k — —£, 



and 



-i — , for k 



!,...,£ + !. 



B. Useful bracket relations 
From the identity 



5 &n,m\ — 3kn&j,m 3jm&ri 



(10) 



we get the relations [E jtk ,E k;n ] = Ej >n , [F j<k ,F kin ] 
-Ej >n , and [E jtk , F ki7l ] = F jt „ and 



[ E j,k,Fj,k] - 2Dj,fc. 



(11) 



The relations above can be interpreted following a "triangle 
rule": the bracket between an operator coupling the states 
Yf 1 and Y£ and an operator coupling the states Y™ and 
Yu, couples the states Y£ and Yff . On the other hand, the 
bracket is zero if two operators couple no common states. 
Moreover, 

[A, -B(£,fc),(£+iji)] = 2(£ + l).F(^fc),(M-i,/i), (12a) 
[A, Fy t k),y+i,h)] = -%(£ + l)E(£ tk ) t ( t+ i th ). (12b) 

From ( [Tol l we find also that 

[E(_l im ),{i+l,m)^{i,m'),{l+l,m'—l)] 

if and only if m! = m or m' = m + 1, with 



and 



£,m),(£+l,m)i 



E, 



,m+l 



),(i+l,m)\ — ^(^,m),(^,m+l)- 



C. Controllability result 

We prove the following result, which allows us to apply 
the abstract controllability criterium obtained in the previsous 
section. We obtain then Theorem 11.11 as a corollary of 
Theorem 13.11 Notice that the conclusions of Theorem 13.11 
allow us to claim more than the required approximately 
controllability, since simultaneous controllability is obtained 
as well. 

Proposition 4.1: The Lie algebra L generated by 
A, Bi,B 2 , B 3 is the whole algebra su(4^ + 4). 

Thanks to the matrix relations obtained in Section IIV-B1 
the proof of the proposition can be easily reduced to the 
proof of the following lemma. 

Lemma 4.2: The Lie algebra L contains the elementary 
matrices 



E, 



(e, k ),(e+i, k +j) 



for k = 



..,£, j = -1,0,1. 



Proof of Lemma 14.21 First, we want to prove that 

{E(t,-j),(t+i,-j) + E(t,j),(t+i,j) I J = 0, • . • C i. (13) 
We use the fact that 

£ 

m——£ 

Indeed, for j = 

m 

[B 3 ,A] = 22 Pt,m[F(£ tm -) f (t + i jm ),A] 
t=—m 

m 

= 2(£ + 1) 2J Pt,mE{t,m),{t+l,m) , 
t=-m 

and, by induction, for j > 1, 

a d 2 J+ 1 A=[B 3) [B 3 ^d%- 1 A}} 

i 

^ Pl*m [ B 3, [Bs, E(i,m),(t+l,m)]] 
m=—£ 

= (-1)^+1)2^ £ pf-l 

m=-£ 

[ B 3, [ P^,h-P(f,fc),(i+l,fe))-E(i,m),(€+l,Tn)]] 

= (-l) J '- 1 (i+l)2 !W - 1 

^ ffyf m [B3, —2pi,mD(t,m),(t+l,m,)] 
rn=~l 

= {-\y{t + i)2 2j 

^ Ppn [ ^ ^,/»i J (i,fe),(«+l,h)i-D(^m),(£+l,m)] 

= (-l)^(^+l)2 2 ^ 1 C%m),(W iro )- 
m— — £ 

Then ( TT3l follows from the fact that p^. m 7^ p^.„ for every 
n m, —to. 
Now note that 

£ 

B 2 — [A,Bx]/(2(£ + 1)) = 2 2J Se,-m^,m),(^+l,fn+l) 
and 

f 

,m),(^+l,m— 1) ■ 

rn— — £ 

Moreover 

[[ ®»m^,m),(^+l,m-l)!^,0),(^+l,0)]>-S(^,0),(^+l,0)] = 

= — &,l[-S(*,0),(*,l)>-^(*>0),(*+l,0)] 

- 3^,o [-5(^+1, -i),(«+i,o))-B^,o),(f+i,o)] 



and, for < fc < £ 

[E ^,-^^,-^,^+1,-^-1)+ 

■ • • + %,-fc+i-B(^,-fe+i),(^+i,-fc) + <&,k-E(*,fe),(*+x,fc-i) + 

+ ■ • ■ + ^(£,-fc),(*+X,-fc) + 

■®(^,]fe),(<+i.*)]'^'(-<.-fc),(-«+i 1 -fc) + ^,fe),(^+i,fe)] 
= —Qi,-k+l [E(£,-k),(£,-k+l) j ^(i,-k),(<+l,-fe)J 

— qt,k[E(t+i,k-i),(e+i,k)i %e,k),(.M-i,k)] 

= <?f,-fc+l-E(f,-fc+l),(f+l,-fc) + <ft,kJ%,fc),(*+X,*-X)- 

Then we get £^ ,-*),(M-x + 
•••> ^,o),(i+x,-i) + -E(&x),(M-x,o) £ L. 
Similarly we can prove that the Lie algebra L contains 

Now, since -E(^, m ),(^+i, m -i) G L and using the relation 

^%,m),(«+i,*»-i)'^(*' m )>(*+ 1 ." 1 ) +-E(^,-m),(«+X,-m) = 
[■^(i+X,77i-X),(i+l,m)i -^(^,m),(^+X,m-X)] = — -^(^,m),(^+X,-m) 

we obtain that £(f, m ),(«+i, TO ) and %,_ m ),(^+i,- m ) belong 
to L for every m = —£, . . . , — 1 

Similarly, -E(f, m ),(^+x,ro) G L implies that 
E(t,m+i),(e+i,m) and ^,-m),(f+x,-m-x) belong to L 
for every to = — £, . . . , — 1 □ 
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